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, Abstract. We construct and analyze a family of coherent states built on sequences 

' of integers originating from the solution of the boson normal ordering problem. These 

sequences generalize the conventional combinatorial Bell numbers and are shown 
to be moments of positive functions. Consequently, the resulting coherent states 
CO ' automatically satisfy the resolution of unity condition. In addition they display such 

\ non-classical fluctuation properties as super-Poissonian statistics and squeezing. 
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Since their introduction in quantum optics many generalizations of 



qh| standard coherent states (CS) have been proposed. 

The main purpose of such generalizations is to account for a full de- 
scription of interacting quantum systems. The conventional CS provide 
^ ' a correct description of a typical non-interacting system, the harmonic 

■ oscillator. One formal approach to this problem is to redefine the stan- 

dard boson creation and annihilation a operators, satisfying [a, a^] = 
1, to A = af{a^a), where the function /(n), n = a^a, is chosen to ad- 
equately describe the interacting problem. Any deviation of f{x) from 
f{x) = const describes a non-linearity in the system. This amounts to 
introducing the modified (deformed) commutation relations ^ [21 El 

[AAt] = [n + l]-[n], (1) 

where the "box" function [n] is defined as [n] = nf'^{n) > 0. Such a way 
of generalizing the boson commutator naturally leads to generalized 
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"nonlinear" CS in the form ([n]! = [0][1] . . . [n], [0] = 1) 

OO yi 

\z)=M-'/\\z\^)Y.^\n), (2) 



71=0 V L J 



which are eigenstates of the "deformed" boson annihilation operator A 

A\z) = z\z). (3) 

The |n)'s are normalized kets in a Fock space. 

It is worth pointing out that such nonlinear CS have been success- 
fully applied to a large class of physical problems in quantum optics 
[^[S]. The comprehensive treatment of CS of the form of Eq. @ can 
be found in [Hill El- 

An essential ingredient in the definition of CS is the completeness 
property (or the resolution of unity condition) [HI E] • A guideline 
for the construction of CS in general has been put forward in as a 
minimal set of conditions. Apart from the conditions of normalizability 
and continuity in the complex label this set reduces to satisfaction 
of the resolution of unity condition. This implies the existence of a 
positive function satisfying 

OO 

cPz\z)W{\z\^){z\= I =Y.\n){n\, (4) 

which reflects the completeness of the set {1^)}. 

In Eq.@ / is the unit operator and \n) is a complete set of or- 
thonormal eigenvectors. In a general approach one chooses strictly 
positive parameters n = 0,1,... such that the state \z) which 

is normalized, {z\z) = 1, is given by 

OO ji 



with normalization 



|2|2" 



which in this note we assume to be a convergent series in \z\'^ for all 
z G C. In view of Eqs. and this corresponds to p{n) = [n]\ or 
[n\ = p{n)/ p{n — 1) for n = 1, 2, . . .. 

Condition (j^ can be shown to be equivalent to the following infinite 
set of equations [7j: 



OO 



W{x) 



dx = p{n), n = 0, 1,..., (7) 
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a Stieltjes moment problem for W{x) = t^jj^- 

Recently considerable progress was made in finding explicit solutions 
of Eq.Q for a large set of p(n)'s, generalizing the conventional choice 
\z)c for which pc{n) = n\ with7V;(x) = (see Refs. [711111111011111112] 
and references therein), thereby extending the known families of CS. 
This progres was faciliated by the observation that when the mo- 
ments form certain combinatorial sequences a solution of the associated 
Stieltjes moment problem may be obtained explicitly {151 . 

In this work we make contact with the combinatorial sequences ap- 
pearing in the solution of the boson normal ordering problem |15l I16j . 
These sequences have the very desirable property of being moments 
of Stieltjes-type measures and so automatically fulfill the resolution of 
unity requirement. It is therefore natural to use these sequences for the 
CS construction, thereby providing a link between the quantum states 
and the combinatorial structures. 

The normal ordering problem for canonical bosons [a, a^] = 1 is re- 
lated to certain combinatorial numbers S{n, k) called Stirling numbers 
of the second kind through |14j 

n 

(flta)" = ^S(n,A;)(at)'=a^ (8) 
fe=i 

with corresponding numbers B{n) = J2k=i 'S'(f^) k) called Bell numbers. 

For integers n, r, s > we define generalized Stirling numbers of the 
second kind Sr,s{n,k) through (r > s): 

ns 

[{a^YaT = ia^r^'-'^ E ^rA^^^ k){a)fa\ (9) 

k=s 

as well as generalized Bell numbers Br^s{n) 

ns 

Br,s{n) = Y.SrAn,k). (10) 

k=s 

For both Sr,si'n,k) and Br^sin) exact and explicit formulas have been 
found Uni EE]. 

In this study we shall only be interested in the subset Br^i{n), for 
which a convenient infinite series representation may be given |15[ I16j: 

This is a generalization of the celebrated Dobihski relation 

fe=0 
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Figure 1. The weight functions Wr,i{x), {x — \z\^), in the resolution of unity for 
r = 2,3,4 (continuous curves) and for r — 1, a Dirac's comb (in the inset), as a 
function of x. 



for conventional Bell numbers -Bi.i(n), see |18j.|19j. For reference we 
quote = 1,2,5,15,52,203,!... 

From Eqs.Q and (|in|) one sees that Bj.^i{n) are integers. Eq. (fTT|l 
gives explicitly B2,i{n) = 1, 3, 13, 73, 501, 4051 . . . ; 53,i(n) = 1,4,25,211,2236,28471,, 
Bi^i{n) = 1, 5, 41, 465, 6721, 117941 . . . etc. Through transformations of 
Ea. (fTT|) one finds that Br^i{n) for r > 1 can be expressed as special 
values of hypergeometric functions of type iF^-i: 

52,i(n) = - iFi(n + l;2;l), (13) 

e 

^ , . 2"-i /2r(n + i) ^ ^ 113 1, 
^^■^("^ = ^H7^^^^^"+2^2'2U) 



etc. 



+n! iF2(n+l;|2;i)) , (14) 



It is essential for our purposes to observe that the integer Br^i{n+1), 
n = 0, 1, . . . is the n-th moment of a positive function Wr^i{x) on the 
positive half-axis. For r = 1 we see from Ea. ()12() that Bi^i{n + 1) is 
the n-th moment of a discrete distribution H^i,i(x) located at positive 
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integers, a so-called Dirac comb: 

POO 

5i,i(n + l) = / X 
Jo 

For every r > 1 a continuous distribution Wr,i{x) will be obtained by 
excising (r — 1)" r(n+^^) from Eq. pij) . performing the inverse Mellin 
transform on it and inserting the result back in the sum of Ea. (|ll|) . 
(Note also that Br^i{0) = r = 2, 3, . . . is no longer integral). In 
this way we obtain 

Bri{n + l)= x'^Wri{x)dx, (16) 
Jo 

which yields for r = 2,3,4: 



1 ^ 6{x 



k) 



1)1 



dx. 



(15) 



a\i(x)=e-'-^-/ih(2^). (17) 
= l/|e-t- (A „^,(1,||) + _L „^,(|,2;|))(,8) 

"'"<^) = iii(iy'''"(^*^^<^'^*»^'ill^)+ 

3'-r.. „f:,(-,-,-;-)+,r(-).„F,(j,-,2;-)), 

In Eqs.((I7I),((IHl),(|ini) Iy{y) and oip(- ••;?/) are modified Bessel and 
hypergeometric functions, respectively. Other Wr,i{x) > for r > 4 
can be generated by essentially the same procedure. 

In Figure Ewe display the weight functions Wr,i{x) for r = 1 ... 4; 
all of them are normalized to one. In the inset the height of the vertical 
line at X = symbolizes the strength of the delta function 5{x — k), 
see Eq.p5|). For further properties of Wi^i{x) and more generally of 
Wr^rix) associated with Eq. (fT(H) . see fTT] . 

A comparison of Eqs.©,® and (|16)1 indicates that the normalized 
states defined through p{n) = Br^i{n + 1) as 

OO 72 

\z)r=M~'/\\z\') Y: , ' In), (20) 

n=0 jBr,l{n + 1) 

with normalization 

OO n 
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Figure 2. Mandel parameters Qr{x), for r = 1 . . . 4, as a function of a; = see 



automatically satisfy the resolution of unity condition of Eq.Q, since 

Nr(x) 



fov Wr,l{x) = TT^^- 



cfz \z)rWr,li\z\'^)r{z\ = / = ^ |"-)("'|- 



(22) 



n=0 



Note that Eg. 1)20(1 is equivalent to Eq.© with the definition [n]r = 
Br,i{n + l)/Br,iin), n = 0,1,2,.... 

Having satisfied the completeness condition with the functions Wr^i{x), 
r = 1, 2, ... we now proceed to examine the quantum-optical fluctuation 
properties of the states \z)r- Prom now on we consider the |n)'s to be 
eigenfunctions of the boson number operator A*" = a^a, i.e. N\n 
The Mandel parameter [7j 

Qr{x) = 



n\n) 



Mr"( 



Mr'{x) 
Mr{x) 



(23) 



allows one to distinguish between the sub-Poissonian (antibunching ef- 
fect, Qr < 0) and super-Poissonian (bunching effect, Qr > 0) statistics 
of the beam. In Pigure[21we display Qr{x) for r = 1 ... 4. It can be seen 
that all the states \z)t in question are super-Poissonian in nature, with 
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Figure 3. The squeezing parameters of Eas. 12411 and 1251 for the coordinate Q (three 
upper curves) and for the momentum P (three lower curves) for different r, as a 
function of Re{z), for r = 1, 2, 3. 



the deviation from Qr = 0, which characterises the conventional CS, 
diminishing for r increasing. 

In Figure we show the behavior of 



and 



5p,,.)..M^. (25) 



which are the measures of squeezing in the coordinate and momentum 
quadratures respectively. In the display we have chosen the section 
along Re{z). All the states \z)r are squeezed in the momentum P 
and dilated in the coordinate Q. The degree of squeezing and dilation 
diminishes with increasing r. By introducing the imaginary part in z 
the curves of Sq{z) and Sp{z) smoothly transform into one another, 
with the identification Sqiia) = Sp{a) and Sp{ia) = SQ{a) for any 
positive a. 

In Figure 0] we show the signal-to-quantum noise ratio j2Uj relative 
to 4[c(2;| A^|z)c] = 4|zp, its value in conventional coherent states; i.e the 
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quantity ar = (Jr — 'i[c{z\N\z)c\, where 

[r{z\Q\z)r? 

with (AQ)^ = r{z\Q'^\r)r — [r{z\Q\r)r]'^ . Again only the section Re{z) is 
shown. We conclude from Figure^that the states \z)r are more "noisy" 
than the standard CS with pdn) = nl . 
In Figure El we give the metric factors 



UJr{x) 



' Mr'{x) 

X ^ — 

J\fr(x) 



(27) 



which describe the geometrical properties of embedding the surface 
of coherent states in Hilbert space, or equivalently a measure of a 
distortion of the complex plane induced by the CS (71 . Here, as far 
as r is concerned, the state \z)i appears to be most distant from the 
\z)c CS for which Uc = 1. 

The use of sequences Br^i{n) to construct CS is not limited to 
the case exemplified by Ea. H2U|) . In fact, any sequence of the form 
Br,i{n + p), P = 0,1,... will also define a set of CS, as then their 
respective weig ht functions wih be v}\\x) = xP-^Wr,i{x) > 0. Will 
the physical properties of CS defined with Pp{n) = Br^i(n + p) depend 
sensitively on p ? Our first guess was that by varying p the salient 
features of the physical results will not change. It was based on a study 
of p{n) = (n + p)! in [7j as a function of p in a somewhat similar 
situation. However, this was not confirmed by actual calculations (and 
in fact could have been envisaged from the outset as V^^^ {x) looks 
quite different from Wr,\{x)). A case in point is that of p = for 
which qualitative differences from Figure |21(p = 1) appear. In Figure 
we present the Mandel parameter for these states. Whereas for r = 1 
the state is still super-Poissonian, for r = 2, 3, 4 one observes novel 
behaviour, namely a crossover from sub- to super-Poissonian statistics 
for finite values of x. Also, as indicated in Figure [7J we note a cross-over 
between squeezing and dilating behaviour for different r. These curious 
features merit further investigation. 

In conclusion, we have used the combinatorial numbers Br^i{n) to 
construct and analyze families of coherent states which automatially 
satisfy the resolution of unity property. This study will be extended 
to other Br^s{n) arising in the problem of boson normal ordering, since 
such sequences, being moments of positive functions are ideally 

suited for the construction of coherent states. 
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Figure 6. Mandcl parameters Qr{x) for r = 1...4, as a function of a; = , for 
states with p{n) = Br,i{n). 



1.5 




Re(z) 



Figure 7. Squeezing parameters for the states with p{n) — Br,i{n) for r = 1, 2, 3, as 
a function of a; = \zf. The subscripts P and Q refer to momentum and coordinate 
variables respectively. 



CombCohStat.tex; 1/02/2008; 7:13; p. 10 



Combinatorial coherent states via normal ordering of bosons 



11 



Acknowledgements 

We thank A. Horzela, M. Mendez, J. Pitman, C. Quesne and A. Vour- 
das for important discussions. 



References 

1. Solomon A.I.: A characteristic functional for deformed photon phenomenology, 
Phys. Lett. A 196 (1994), 29-34. 

2. Man'ko, V.I., Marmo, G., Sudarshan, E.C.G. and Zaccaria, F.: f-Oscilators and 
nonlinear coherent states, Phys.Scr. 55 (1997), 528-541. 

3. Man'ko, V.I. and Vilela Mendes, R.; On the nonlinearity interpretation of q- 
and /-deformations and some applications, J. Phys. A 31 (1998), 6037-6044. 

4. Matos Filho, R.L. and Vogel, W.: Nonlinear coherent states, Phys. Rev. A 54 
(1996), 4560-4563. 

5. Kis, Z., Vogel, W. and Davidovich, L.: Nonlinear coherent states of trapped- 
atom motion, Phys. Rev. A 64 (2001), 033401-033411. 

6. Klauder, J.R.: Continuous-representation theory, J. Math. Phys. 4 (1963), 
1055-1073. 

7. Klauder, J.R., Penson, K.A., Sixdeniers, J.M.: Constructing coherent states 
through solutions of Stieltjes and Hausdorff moment problems, Phys. Rev. A 
64 (2001), 013817-013835. 

8. Penson, K.A., Solomon, A.I.: Mittag-LefHer coherent states, J. Phys. A 32 
(1999), 7543-7563; Sixdeniers, J.M. and Penson, K.A.: On the completeness 
of coherent states generated by binomial distribution, J. Phys. A 33 (2000), 
2907-2916; Sixdeniers, J.M., Penson, K.A. and Klauder, J.R.: Tricomi coherent 
states. Int. J. Mod. Phys. B 15 (2001), 4231-4243. 

9. Quesne, C: Generalized coherent states associated with the Cx-extended os- 
cillator, Ann. Phys. (N.Y.) 293 (2001), 147-188; Quesne C: New q-deformcd 
coherent states with an explicitly known resolution of unity, J. Phys. A 35 
(2002), 9213-9226. 

10. Vourdas, A. and Wiinsche, A.: Resolutions of the identity in terms of line 
integrals of SU(1, 1) coherent states, J. Phys. A 31 (1998), 9341-9352; Vourdas, 
A.: The growth of Bargmann functions and the completeness of sequences of 
coherent states, J. Phys. A 30 (1997), 4867-4876. 

11. Horowski, M., Odzijewicz, A. and Tereszkiewicz, A.: Some integrable systems 
in nonlinear quantum optics, J. Math. Phys. 44 (2003), 480-506. 

12. Novaes, M. and Gazeau, J. P.: Multidimensional generalized coherent state, J. 
Phys. A 36 (2003), 199-212. 

13. Penson, K.A. and Solomon, A.I.: Coherent states from combinatorial sequences. 
In: E. Kapuscik and A. Horzela (eds). Proceedings of the 2nd International 
Symposium on Quantum Theory and Symmetries, World Scientific, Singapore, 
2001, pp. 527-530. 

14. Katriel, J.: Combinatorial aspects of boson algebra, Lett. Nuovo Cimento 10 
(1974), 565-567. 

15. Blasiak, P., Penson, K.A. and Solomon, A. I.: The general boson normal 
ordering problem, Phys. Lett. A 309 (2003), 198-205. 



CombCohStat.tex; 1/02/2008; 7:13; p. 11 



12 Blasicik P, Penson K A and Solomon A I 

16. Blasiak, P., Penson, K.A. and Solomon, A. I.: The boson normal ordering 
problem and generalized Bell numbers, Ann. Comb. 7 (2003), 127-139. 

17. Blasiak, P., Penson, K.A. and Solomon, A.I.: Dobinski-type relations and the 
log-normal distribution, J. Phys. A 36 (2003), L273-L278. 

18. Wilf, H.S.: Generatingfunctionology, Academic Press, New York, 1994. 

19. Pitman, J.: Some probabilistic aspects of set partitions, Amer. Math. Monthly 
104 (1997), 201-209. 

20. Yuen, H.P.: Two-photon coherent states of the radiational field, Phys. Rev. A 
13 (1976), 2226-2243. 



CombCohStat.tex; 1/02/2008; 7:13; p. 12 



